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At ultra-high density, matter is expected to form a degenerate Fermi gas of quarks 
in which there is a condensate of Cooper pairs of quarks near the Fermi sur- 
face: color superconductivity. In this chapter we review some of the underlying 
physics, and discuss outstanding questions about the phase structure of ultra- 
dense quark matter. We then focus on describing recent results on the crystalline 
color superconducting phase that may be the preferred form of cold, dense but 
not asymptotically dense, three-flavor quark matter. The gap parameter and free 
energy for this phase have recently been evaluated within a Ginzburg-Landau 
approximation for many candidate crystal structures. We describe the two that 
are most favorable. The robustness of these phases results in their being favored 
over wide ranges of density. However, it also implies that the Ginzburg-Landau 
approximation is not quantitatively reliable. We describe qualitative insights into 
what makes a crystal structure favorable which can be used to winnow the pos- 
sibilities. We close with a look ahead at the calculations that remain to be done 
in order to make quantitative contact with observations of compact stars. 



1. Introduction 

The exploration of the phase diagram of matter at ultra- high temperature or density- 
is an area of great interest and activity, both on the experimental and theoretical 
fronts. Heavy-ion colliders such as the SPS at CERN and RHIC at Brookhaven have 
probed the high-temperature region, creating and studying the properties of quark 
matter with very high energy density and very low baryon number density similar to 
the fluid which filled the universe for the first microseconds after the big bang. In this 
paper we discuss a different part of the phase diagram, the low-temperature high- 
density region. Here there are as yet no experimental constraints, and our goal is to 
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understand the properties of matter predicted by QCD well enough to be able to use 
astronomical observations of neutron stars to learn whether these densest objects 
in the current universe contain quark matter in their core. We expect cold, dense, 
matter to exist in phases characterized by Cooper pairing of quarks, i.e. color su- 
perconductivity, driven by the Bardccn-Cooper-Schricffcr (bcs; 1 mechanism. The 
BCS mechanism operates when there is an attractive interaction between fermions 
at a Fermi surface. The QCD quark-quark interaction is strong, and is attractive in 
many channels, so we expect cold dense quark matter to generically exhibit color 
superconductivity. Moreover, quarks, unlike electrons, have color and flavor as well 
as spin degrees of freedom, so many different patterns of pairing are possible. This 
leads us to expect a rich phase structure in matter beyond nuclear density. 

Calculations using a variety of methods agree that at sufficiently high density, 
the favored phase is color-flavor-locked (CFL) color-superconducting quark matted 
(for reviews, see Ref. |3J). However, there i s sti ll uncertainty over the nature of the 
next phase down in density. Previous wor had suggested that when the density 
drops low enough so that the mass of the strange quark can no longer be neglected, 
there is a continuous phase transition from the CFL phase to a new gapless CFL 
(gCFL) phase, which could lead to observable consequences if it occurred in the cores 
of neutron starsP^ However, it now appears that some of the gluons in the gCFL 
phase have i maginary Meissner masses , ind icating an instability towards a lower- 
energy phase! ' x v l LL l lz l i-v l 14 l l \> i l Analysis in the vicinity of the unstable 
gCFL phase cannot determine the nature of the lower-energy phase that resolves 
the instability. However, the instability is telling us that the system can lower its 
energ y by turning on currents, suggesting that the crystalline color superconducting 
phasef^in which the condensate is modulated in space in a way that can be thought 
of as a sum of counterpropagating currents, is a strong candidate. 



2. Review of color superconductivity 
2.1. Color superconductivity 

The fact that QCD is asymptotically free implies that at sufficiently high density 
and low temperature, there is a Fermi surface of weakly-interacting quarks. The 
interaction between these quarks is certainly attractive in some channels (quarks 
bind together to form baryons), so we expect the formation of a condensate of 
Cooper pairs. We can see this by considering the grand canonical potential F = 
E — uN, where E is the total energy of the system, /i is the chemical potential, and 
N is the number of quarks. The Fermi surface is defined by a Fermi energy Ep = a, 
at which the free energy is minimized, so adding or subtracting a single particle 
costs zero free energy. Now switch on a weak attractive interaction. It costs no free 
energy to add a pair of particles (or holes), and if they have the right quantum 
numbers then the attractive interaction between them will lower the free energy of 
the system. Many such pairs will therefore be created in the modes near the Fermi 
surface, and these pairs, being bosonic, will form a condensate. The ground state 
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will be a superposition of states with all numbers of pairs, breaking the fermion 
number symmetry. 

A pair of quarks cannot be a color singlet, so the resulting condensate will break 
the local color symmetry SU(3) co \ OT . The formation of a condensate of Cooper pairs 
of quarks is therefore called "color superconductivity" . The condensate plays the 
same role here as the Higgs condensate does in the standard model: the color- 
superconducting phase can be thought of as the Higgs phase of QCD. 

2.2. Highest density: Color-flavor locking (CFL) 

It is by now well-established that at sufficiently high densities, where the up, down 
and strange quarks can be treated on an equal footing and the disruptive effects 
of the strange quark mass can be neglected, quark matter is in the color-flavor 
locked (CFL) phase, in which quarks of all three colors and all three flavors form 
conventional Cooper pairs with zero total momentum, and all fcrmionic excitations 
are gapped, with th e gap parameter Ao ~ 10-100 MeV.EElThis has been c onfirmed 
by both NJlPEOEland gluon-mediated interaction calculat ions EIE1E1E1 The CFL 
pairing pattern iJ2 

(q?C l5 qP) = A ( K + l)6?6f + A (k - 1)5 f 5? = A e a ^ N e ijN + A k(- ■ ■ ) 

[St/(3)coior] x SU(3) L x SU(3) R xU(1) b -► SU(3) c+ l+r xZ 2 (1) 
3 [U(1)q] D [U{1)q] 

Color indices a,/3 and flavor indices i,j run from 1 to 3, Dirac indices are sup- 
pressed, and C is the Dirac charge-conjugation matrix. The term multiplied by k 
corresponds to pairing in the (6s, 6s), which although not energetical ly favored 
breaks no additional symmetries and so K is in general small but nonzero02E32] 
The Kronecker S's connect color indices with flavor indices, so that the condensate 
is not invariant under color rotations, nor under flavor rotations, but only under 
simultaneous, equal and opposite, color and flavor rotations. Since color is only a 
vector symmetry, this condensate is only invariant under vector flavor+color rota- 
tions, and breaks chiral symmetry. The features of the CFL pattern of condensation 
ar<P 

— The color gauge group is completely broken. All eight gluons become mas- 
sive. This ensures that there are no infrared divergences associated with 
gluon propagators. 

— All the quark modes are gapped. The nine quasiquarks (three colors times 
three flavors) fall into an 8 © 1 of the unbroken global SU(3). Neglecting 
corrections of order k (see l|T)l) the octet and singlet quasiquarks have gap 
parameter Ao and 2Ao respectively. 

— A rotated electromagnetic gauge symmetry (generated by "Q" ) survives un- 
broken. The single remaining massless gauge boson is a linear combination 
of the original photon and one of the gluons. 
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— Two global symmetries are broken, the chiral symmetry and baryon num- 
ber, so there are two gauge-invariant order parameters that distinguish the 
CFL phase from the QGP, and corresponding Goldstone bosons which are 
long-wavelength disturbances of the order parameter. When the light quark 
mass is nonzero it explicitly breaks the chiral symmetry and gives a mass 
to the chiral Goldstone octet, but the CFL phase is still a superffuid, dis- 
tinguished by its spontaneous baryon number breaking. 



2.3. Less dense quark matter: stresses on the CFL phase 

The CFL phase is characterized by pairing between different flavors and different 
colors of quarks. We can easily understand why this is to be expected. Firstly, 
the QCD interaction between two quarks is most attractive in the channel that is 
antisymmetric in color (the 3). Second ly, pairing tends to be stronger in channels 
that do not break rotational symmetry p^^7 l ^8 l ^9Ki Dl so we expect the pairs to 
be spin singlets, i.e. antisymmetric in spin. Finally, fcrmionic antisymmetry of the 
Cooper pair wavefunction then forces the Cooper pair to be antisymmetric in flavor. 

Pairing between different colors/flavors can occur easily when they all have the 
same chemical potentials and Fermi momenta. This is the situation at very high den- 
sity, where the strange quark mass is negligible. However, even at the very center of 
a compact star the quark number chemical potential [i cannot be much larger than 
500 MeV, meaning that the strange quark mass M s (which is density dependent, ly- 
ing somewhere between its vacuum current mass of about 100 MeV and constituent 
mass of about 500 MeV) cannot be neglected. Furthermore, bulk matter, as rele- 
vant for a c ompact star, must be in weak equilibrium and must be electrically and 
color neutra l I I I I I I (possibly via mixing of oppositely charged phases). 
All these factors work to separate the Fermi momenta of the three different flavors 
of quarks, and thus disfavor the cross-species BCS pairing that characterizes the 
CFL phase. If we imagine beginning at asymptotically high densities and reducing 
the density, and suppose that CFL pairing is disrupted by the heaviness of the 
strange quark before color superconducting quark matter is superseded by baryonic 
matter, the CFL phase must be replaced by some phase of quark matter in which 
there is less, and less symmetric, pairing. 

In the next few subsections we give a quick overview of the expected phases of 
real-world quark matter. We restrict our discussion to zero temperature because the 
critical temperatures for most of the phases that we discuss are expected to be of 
order 10 MeV or higher, and the core temperature of a neutron star is believed to 
drop below this value within minutes (if not seconds) of its creation in a supernova. 

2.4. Kaon condensation: the CFL-K phase 

Bedaque and SchafeJ232H showed that when the stress is not too large (high den- 
sity), it may simply modify the CFL pairing pattern by inducing a flavor rotation 
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of the condensate which can be interpreted as a condensate of mesons, i.e. the 

neutral anti-strange Goldstonc bosons associated with the chiral symmetry break- 
ing. This is the "CFL-KO" phase, which breaks isospin. The if°-condensate can 
easily be suppressed by instanton effectspHl but if these are ignored then the kaon 
condensation occurs for M s > m 1 / 3 A 2 / 3 for light (u and d) quarks of mass m. This 
was demonstrated using an effective theory o f the G oldstone bosons, but with some 
effort can also be seen in an NJL calculationEHEIl 



2.5. The gapless CFL phase 

The nature of the next significant transition has been studied in NJL model cal- 
culations which ign ore the ^condensat ion in the CFL phase and which assume 
spatial homogeneity N^2 | 4? | 44 | 4^ | 4( j ? | 47 | nas k een found that the phase struc- 
ture depends on the strength of the pairing. If the pairing is sufficiently strong (so 
that A ~ 100 MeV where A is what the CFL gap would be at /i - 500 MeV 
if M s were zero) then the CFL phase survives all the way down to the transition 
to nuclear matter. For a wide range of parameter values, however, we find some- 
thing more interesting. We can make a rough quantitative analysis by expanding in 
powers of M 2 /fi 2 and A///, and ignorin g th e fact that the effective strange quark 
mass may be different in different phases!""' Such an analysis shows that, within a 
spatially homogeneous ansatz, as we come down in density we find a transi tion at 
/i ps ±M%/A from CFL to another phase, the gapless CFL phase (gCFL).EEIln 
this phase, quarks of all three colors and all three flavors still form ordinary Cooper 
pairs, with each pair having zero total momentum, but there are regions of momen- 
tum space in which certain quarks do not succeed in pairing, and these regions are 
bounded by momenta at which certain fermionic quasiparticles are gapless. This 
variation on BCS pairing — in which the same species of fermions tha t pair f eature 
gapless quasiparticles — was first proposed for two flavor quark matteJ^ESI an d i n 
an atomic physics context!^ 

For Ml) ' n > 2Ao, the CFL phase has higher free energy than the gCFL phase. 
This follows from the energetic balance between the cost of keeping the Fermi sur- 
faces together and the benefit of the pairing that can then occur. The leading effect 
of M s is like a shift in the chemical potential of the strange quarks, so the bd and gs 
quarks feel "effective chemical potentials" /uff = // — §A*8 an d Mgf = A* + — "^7- 
In the CFL phase ^8 = -M s 2 /(2/i)ESEi S0 ^ff _ ^ = M*/fi. The CFL phase will 
be stable as long as the pairing makes it energetically favorable to maintain e quali ty 
of the bd and gs Fermi momenta, despite their differing chemical potentials! I I It 
becomes unstable when the energy gained from turning a gs quark near the common 
Fermi momentum into a bd quark (namely / /x) exceeds the cost in lost pairing 
energy 2Arj. So the CFL phase is stable when 

M 2 

— ^ < 2A . (2) 



February 2, 2008 6:19 WSPC/Trim Size: 9.75in x 6.5in for Review Volume beyond'BCS 



6 




Fig. 1. Gap parameter A versus M s 2 //i for the CFL and gCFL phases as well as for crystalline 
color superconducting phases with two different crystal structures described in Scction l3.6l These 
gap parameters have all been evaluated in a NJL model with the interaction strength chosen such 
that the CFL gap parameter is 25 MeV at /^t = 0. In the gCFL phase, A ud > A us > A ds . The 
calculation of the gap parameters in the crystalline phases has been done in a Ginzburg-Landau 
approximation. Recall that the splitting between Fermi surfaces is proportional to M% //^, and that 
small (large) / fj, corresponds to high (low) density. 



Including the effects of A'°-condensation expa nds the r ange of Mj* j [i in which the 
CFL phase is stable by a few tens of percent! I I I For larger M^/fi, the CFL 
phase is replaced by some new phase with unpaired bd quarks, which cannot be 
neutral unpaired or 2SC quark matter because the new phase and the CFL phase 
must have the same free energy at the critical M^/fj, = 2Ar> 

The obvious approach to finding this phase is to perform a NJL model calculation 
with a general ansatz for the pairing that includes differences between the flavors, 
for example by allowing different pairing strengths A uc i, Ad s , A us . This was done 
in Refs. 0] [SJ and the resultant "gCFL" phase was described in detail. In Figs. ^ 
and [2] we show the results. The gCFL phase takes over from CFL at / [i ss 2Ao, 
and remains favored beyond the value M^/fi ss 4Ao at which the CFL phase would 
become unfavored. 



2.6. Beyond gapless CFL 

The arguments above led us to the conclusion that the favored phase of quark 
matter at the highest densities is the CFL phase, and that as the density is de- 
creased there is a transition to another color superconducting phase, the gapless 
CFL phase. However, it turns out that the gCFL phase is itself unstable, and that 
in the regime where the gCFL phase has a lower free energy than both the CFL 
phase and unpaired quark matter, there must be some other phase with lower free 
energy still. The nature of that phase has to date not been determined rigorously, 
and in subsequent subsections we shall discuss various possibilities. However, the 
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Fig. 2. Free energy Q relative to that of neutral unpaired quark matter versus / fi for the CFL, 
gCFL and crystalline phases whose gap parameters are plotted in Fig. ^ Recall that the gCFL 
phase is known to be unstable, meaning that in the regime where the gCFL phase free energy is 
plotted, the true ground state of three-flavor quark matter must be some phase whose free energy 
lies below the dashed line. We see that the three-flavor crystalline color superconducting quark 
matter phases with the most favorable crystal structures that we have found, namely 2Cube45z 
and CubeX described in 1311 and 1331 . have sufficiently robust condensation energy (sufficiently 
negative Q) that they are candidates to be the ground state of matter over a wide swath of / fi, 
meaning over a wide range of densities. 



crystalline color superconducting phases analyzed very recently in Rcf. whose 
free energies are shown in Fig.[21appcar particularly robust, offering the possibility 
that over a wide range of densities they are the ground state of dense matter. We 
shall therefore discuss these phases at much greater length in Section |21 

In all the contexts in which they have been investigated (whether in atomic 
physics or the g2SC phase in two- flavor quark matter or the gCFL phase) the gapless 
paired state turns out to suffer from a "magnetic i nstability" - it can lower its ener gy 
by the formation of counter-propagating currents! ********* * ' l 10 l l n 
the atomic physics context, the resolution of the instability is phase separation 
into macroscopic regions of two phases in one of which sta ndard BCS pairing oc- 
curs and in the other of which no pairing occurs P S^^ I ^ I Phase separation into 
electrically charged but color neutral phases is also a possibility in two-flavor quark 
matter.^l In three-flavor quark matter, where the instability of the gCFL phase has 
been established in Refs.03El phase coexistence would require coexisting compo- 
nents with opposite color charges, in addition to opposite electric charges, making 
it very unli kely t hat a phase separated solution can have lower energy than the 
gCFL phase! * I Furthermore, color superconducting phases which are less sym- 
metric than the CFL phase but still involve only conventional BCS pairing, for 
example th e much-st udied 2SC phase in which only two colors of up and down 
quarks pair | 57 | 58 | 59 | includin g also m any other possibilities ESI cannot be the 
resolution of the gCFL instability! '"^'^ It seems likely, therefore, that a ground 
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Fig. 3. Dispersion relations of the lightest quasiquark excitations in the gCFL phase, at /i = 
500 McV, with M a = 200 MeV and interaction strength such that the CFL gap parameter at M s = 
would be Ao = 25 MeV. Note that in there is a gapless mode with a quadratic dispersion relation 
(energy reaching zero at momentum p\ u 2 ) as well as two gapless modes with more conventional 
linear dispersion relations. 



state with counter-propagating currents is required. This could take the form of 
a crystalline color supcrconducto J 19 ' ^ 1 ' ^ 2 ' 6 ^ ' ^^^^^ ' ^^^^^^?'' ! ?! ! ^^ — the QCD 

analogue of a fo rm of non-BCS pairing first considered by Larkin, Ovchinnikov, 
Fulde and Fe rrellPlpr, given that the CFL phase itself is likely augmented by kaon 
condensation I2Z1 38 | 4() | 41 | ^ cou id take the form of a phase in which a CFL kaon con- 
densate carries a current in one direction balanced by a counter-p ropaga ting current 
in the opposite direction carried by gapless quark quasipart icles ! ' 

The instability of the gCFL phase appears to be related to one of its most 
interesting features, namely the presence of gapless fermionic excitations around 
the ground state. These are illustrated in Fig. |3 which shows that there is one 
mode (the bu-rs quasiparticle) with an unusual quadratic dispersion relation, which 
is expected to give rise to a parametrically enhanced heat capacity and neutrino 
emissivity and anomalous transport properties.^ The instability manifests itself in 
imaginary Meissner masses Mm for some of the gluons. M M is the low-momentum 
current-current two-point function, and M M /(g 2 A 2 ), with g the gauge coupling, 
is the coefficient of the gradient term in the effective theory of small fluctuations 
around the ground-state condensate. The fact that we find a negative value when 
the quasiparticles are gapless indicates an instability towards spontaneo us b reak- 
ing of translational invariance. Calculations in a simple two-species mode^^ show 
that imaginary Mm is generically associated with the presence of gapless charged 
fermionic modes. 

The earliest calculations for the three-flavor case show that even a very simple 
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ansatz for the crystal structure yields a crystalline color superconducting state that 
has lo wer free energy than gCFL in the region where the gCFL^unpaired transition 
occurs ESllZll it is reasonable, based on what was found in the two-flavor case in 
Ref. EH1 to expect that when the full space of crystal structures is explored, the 
crystalline color superconducting state will be preferred to gCFL over a much wider 
range of the stress parameter M 2 '/( Ap/i), and this expectation has recently been 
confirmed by explicit calculation!^! It is conceivable that the whole gCFL region 
is actually a crystalline region, but the results shown in Fig. [21 that we describe in 
Section [21 suggest that there is still room for other possibilities (like the current- 
carrying meson condensate) at the highest densities in the gCFL regime. 

An alternative explanation of the consequences of the gCFL instability was 
advanced by HonJ23 (see also Ref. I49f) : since the instability is generically associated 
with the presence of gapless fermionic modes, and the BCS mechanism implies that 
any gapless fermionic mode is unstable to Cooper pairing in the most attractive 
channel, one might expect that the instability will simply be resolved by "secondary 
pairing" . This means the formation of a (qq) condensate where q is either one of 
the gapless quasiparticles whose dispersion relation is shown in Fig. [21 After the 
formation of such a secondary condensate, the linear gapless dispersion relations 
would be modified by "rounding out" of the corner where the energy falls to zero, 
leaving a secondary energy gap A s , which renders the mode gapped, and removes the 
instability. In the case of the quadratically gapless mode there is a greatly increased 
density of states at low energy (in fact, the density of states diverges as E^ 1 / 2 ), so 
Hong calculated that the secondary pairing should be much stronger than would 
be predicted by BCS theory, and he specifically predicted A s cx G 2 for coupling- 
strength G s , as compared with the standard BCS result A cx exp(— const/G s ). 

This possibility was worked out in a two-species NJL model in Ref. [2] This 
allowed a detailed exploration of the strength of secondary pairing. The calculation 
confirmed Hong's prediction that in typical secondary channels A s cx G 2 S . However, 
in all the secondary channels that were analyzed it was found that the secondary 
gap, even with this enhancement, is from ten to hundreds of times smaller than 
the primary gap at reasonable values of the secondary coupling. This shows that 
that secondary pairing does not generically resolve the magnetic instability of the 
gapless phase, since it indicates that there is a temperature range A s C T « A p 
in which there is primary pairing (of strength A p ) but no secondary pairing, and 
at those temperatures the instability problem would arise again. 

2.7. Crystalline pairing 

The pairing patterns discussed so far have been translationally invariant. But in the 
region of parameter space where cross-species pairing is excluded by stresses that 
pull apart the Fermi surfaces, one expects a position-d ependent pairin g known as 
the crystalline color superconducting phase! 72 l 19 l 61 l 62 l 63 l 64 l 65 l 66 i 67 l 68 i 69 l 7t) l 71 l ; ^ l 

This arises because one way to achieve pairing between different flavors while acco- 
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Fig. 4. Pictorial representation of CSL pairing, 2SC+CSL pairing, and 2SC+1SC pairing in 
neutral quark matter. See Ref. 1751 



modating the tendency for the Fermi momenta to separate is to allow Cooper pairs 
with nonzero total momentum yielding pairing over parts of the Fermi surfaces. We 
describe this phase extensively in Section [3J As we will discuss there, this phase 
may resolve the g'CFL phase's stability problems. 



2.8. Single-flavor pairing 

At densities that are so low that M s puts such a significant stress on the pairing 
pattern that even the crystalline phase becomes unfavored, wc expect a transition 
to a phase with no cross-species pairing at all. (This regime will only arise if Ao 
is so small that very large values of Af^/(/iAo) can arise without /1 being taken so 
small that nuclear matter becomes favored.) 

"Unpaired" quark matter. In most NJL studies, matter with no cross-species pair- 
ing at all is described as as "unpaired" quark matter. However, it is well known 
that there are attractive channels for a single flavor pairing, although they are 
much weaker than the 2SC and CFL channels! I | ^o | zt) | oU| Calculations using 
NJL models and single-gluon exch ang e agree that the favored phase in this case is 
the color-spin-locked (CSL) phasej^ in which all three colors of each flavor, with 
each pair of colors correlated with a particular direction for the spin. This phase 
does not break rotational symmetry. See first panel of Fig. 0] 
"Unpaired" species in 2SC quark matter. 

If the stress is large enough to destroy CFL pairing, and furthermore is large 
enough to preclude cross-species pairing in which the Cooper pairs have nonzero 
momentum and form a crystalline color superconducting phase, we can ask whether 
there is any pattern of BCS pairing which can persist in neutral quark matter 
with arbitrarily large splitting between Fermi surfaces. The answer to this is no, as 
demonstrated via an exhaustive search in Ref. I6UI Although unlikely, it remains a 
possibility that up and down quarks with two colors could pair in the 2SC pattern. 
This does not happen in NJL models in which M s is a parameter. However, in 
models in which one instead solves self-consistently for the quark masses, the 2SC 
phase tends to become more robust. It remains to be seen whether this relative 
strengthening of the 2SC phase in such models is sufficient to allow it to exist at 
values of / fi which are larger than those where the crystalline phases shown in 
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Fig. |21 and discussed at length in Section are favored. 

If there is a regime in which the 2SC phase survives, this leaves the blue quarks 
unpaired. In that case one might expect a "2SC+CSL" pattern, illustrated in the 
second panel of Fig. 01 which would again be rotationally symmetric. However, the 
2SC pattern breaks the color symmetry, and in order to maintain color neutrality, a 
color chemical potential is generated, which also affects the unpaired strange quarks, 
splitting the Fermi momentum of the blue strange quarks away from that of the 
red and green strange quarks. This is a small effect, but so is the CSL pairing gap, 
and NJL model calculations indicate that the color chemical potential typically 
destroys CSL pairing of the strange quarksl^l The system falls back on the next 
best alternative, which is spin-I pairing of the red and green strange quarks (third 
panel of Fig. 0}. 

Single-flavor pairing phases have much lower critical temperatures than multi- 
flavor phases like the CFL or crystalline p hases, perhaps as large as a few MeV, 
more typically in the eV to many keV range EEEZEHEHEDl so they are expected to 
play a role late in the life of a neutron star. 

2.9. Mixed Phases 

Another way for a system to deal with a stress on its pairing pattern is phase sepa- 
ration. In the context of quark matter this corresponds to relaxing the requirement 
of local charge neutrality, and requiring neutrality only over long distances, so we 
allow a mixture of a positively charged and a negatively charged phase, with a 
common pressure and a common value of the electron chemical potential \i e that 
is not equal to the neutrality value for either phase. Such a mixture of nuclear and 
CFL quark matter was studied in Ref. [THJ In three-flavor quark matter it has been 
found that as long as we require local color neutrality such mixed phases are not 
the favored response to the stress imposed by the st ran ge quark mass!^M2l Phases 
involving color charge separation have been studied^ but it seems likely that the 
energy cost of the color-electric fields will disfavor them. 

3. The crystallography of three-flavor quark matter 

In this section we review recent work by Rajagopal and Sharma!^ Here, we focus 
on explaining the context, the model, the results and the implications. We do not 
describe the calculations. 



3.1. Introduction and context 

The investigation of crystalline color superconductivity in three-flavor QCD was ini- 
tiated in Ref. EH Although such phases seem to be free from magnetic instability!^ 
it remains to be determined whether such a phase can have a lower free energy 
than the meson current phase, making it a possible res olution to the gCFL in- 
stability. The simplest "crystal" structures do not suffice EMLLI b u t experience in 
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the two-flavor context!*^! suggests that realistic crystal structures constructed from 
more plane waves will prove to be qualitatively more robust. The results of Ref. 1521 
confirm this expectation. 

Determining the favored crystal structure(s) in the crystalline color supercon- 
ducting phase(s) of three-flavor QCD requires determining the gaps and comparing 
the free energies for very many candidate structures, as there are eve n m ore possi- 
bilities than the many that were investigated in the two-flavor context!^ As there, 
we shall make a Ginzburg-Landau approximation. This approximation is controlled 
if A Ao, where A is the gap parameter of the crystalline color superconducting 
phase itself and Ao is the gap parameter in the CFL phase that would occur if M s 
were zero. We shall find that the most favored crystal structures can have A/Arj as 
large as ~ 1/2, meaning that we are pushing the approximation hard and so should 
not trust it quantitatively Earlier work on a particularl y sim ple one parameter fam- 
ily of "crystal" structures in three-flavor quark matterf^ simple enough that the 
analysis could be done both with and without the Ginzburg-Landau approximation, 
demonstrates that the approximation works when it should and that, at least for 
very simple crystal structures, when it breaks down it always underestimates the 
gap A and the condensation energy. Furthermore, the Ginzburg-Landau approxima- 
tion correctly determines which crystal structure among the one parameter family 
that we analyzed in Ref. [2] has the largest gap and lowest free energy. 

The underlying microscopic theory that we use is an NJL model in which the 
QCD interaction between quarks is replaced by a point-like four-quark interaction 
with the quantum numbers of singlc-gluon exchange, analyzed in mean field the- 
ory. This is not a controlled approximation. However, it suffices for our purposes: 
because this model has attraction in the same channels as in QCD, its high density 
phase is the CFL phase; and, the Fermi surface splitting effects whose qualitative 
consequences we wish to study can all be built into the model. Note that we shall 
assume throughout that Ao <C \i. This weak coupling assumption means that the 
pairing is dominated by modes near the Fermi surfaces. Quantitatively, this means 
that results for the gaps and condensation energies of candidate crystalline phases 
are independent of the cutoff in the NJL model when expressed in terms of the CFL 
gap Ao: if the cutoff is changed with the NJL coupling constant adjusted so that Ao 
stays fixed, the gaps and condensation energies for the candidate crystalline phases 
also stay fixed. This makes the NJL model valuable for making the comparisons 
that are our goal. 

We shall consider crystal structures in which there are two condensates 



As in Refs. HSU 13 we neglect ids) pairing because the d and s Fermi surfaces are 
twice as far apart from each other as each is from the intervening u Fermi surface. 



(ud) ~ A 3 

exp (2zq3 • r) 



a 




(3) 



a 
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Were we to set A 2 to zero, treating only (ud) pairing, we would recover the two- 
flavor Ginzburg-Landau analysis of Ref.EHl There, it was found that the best choice 
of crystal structure was one in which pairing occurs for a set of eight q 3 's pointing at 
the corners of a cube in momentum space, yielding a condensate with face-centered 
cubic symmetry. The analyses of three-flavor crystalline color superconductivity in 
Refs. 1711 introduced nonzero A 2 , but made the simplifying ansatz that pairing 
occurs only for a single q 3 and a single q 2 . We consider crystal structures with up 
to eight qg's and up to eight q 2 's. 

We evaluate the free energy f2(A 2 ,A3) for each three-flavor crystal structure 
in a Ginzburg-Landau expansion in powers of the A's. We work to order A 2 A| 
with p + q = 6. At sextic order, we find that fl(A, A) is positive for large A for 
all the crystal structures that we investigate!^ This is in marked contrast to the 
result that many two-flavor crystal structures have negative sextic terms, with free 
energies that are unbounded from below when the Ginzburg-Landau expansion is 
stopped at sextic orderP^l Because we find positive sextic terms, we are able to 
use our sextic Ginzburg-Landau expansion to evaluate A and f2(A,A) for all the 
structures that we analyze. 

We can then identify the two most favorable crystal structures, 2Cube45z and 
CubeX, which are described further in Section 13.61 To a large degree, our argu- 
ment that these two structures are the most favorable relies only on two quali- 
tative factors, which we explain in Section 13.51 CubeX and 2Cube45z have gap 
parameters that are large: a third to a half of Ao. We have already discussed the 
implications of this for the Ginzburg-Landau approximation in Section 13.11 above. 
The large gaps naturally lead to large condensation energies, easily 1/3 to 1/2 
of that in the CFL phase with M s = 0, which is 3Aq^ 2 /it 2 . This is remarkable, 
given the only quarks that pair are those lying on (admittedly many) rings on the 
Fermi surfaces, whereas in the CFL phase with M s = pairing occurs over the 
entire u, d and s Fermi surfaces. The gapless CFL (gCFL) phase provides a use- 
ful comparison at nonzero M s . For 2Ao < M 2 /fi < 5. 2An. m odel analyses that 
are restricted to isotropic phases predict a gCFL phase,™"^ finding this phase 
to have lower free energy than either the CFL phase or unpaired quark matter, 
as shown in Fig. However, this phase is unstabl e to the formation of current- 
carrying condensates ! ' ' ' ^ ' ' I i4 l 1 v l 1 9 1 1 < 1 1 v l and so it cannot be the ground 
state. The true ground state must have lower free energy than that of the gCFL 
phase, and for this reason the gCFL free energy provides a useful benchmark. We 
find that three-flavor crystalline color superconducting quark matter with one or 
other of the two crystal structures that we argue are most favorable has lower free 
energy (greater condensation energy) than CFL quark matter, gCFL quark matter, 
and unpaired quark matter for 

2.9A < — < 10.4A . (4) 
A* 

(See Fig. [2] above.) This window in parameter space is in no sense narrow. Our 
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results therefore indicate that three-flavor crystalline quark matter will occur over 
a wide range of densities. (That is unless the pairing between quarks is so strong, i.e. 
Ao so large making Ml / Ao so small, that quark matter is in the CFL phase all the 
way down to the density at which quark matter is superseded by nuclear matter.) 
However, our results also indicate that unless the Ginzburg-Landau approximation 
is underestimating the condensation energy of the crystalline phase by about a factor 
of two, there is a fraction of the "gCFL window" (with 2A < M 2 /^ < 2.9A , in 
the Ginzburg-Landau approximation) in which no crystalline phase has lower free 
energy than the gCFL phase. This is thus the most likely regime in which to find 
the current-carrying meson condensates of Refs. El DSI 

3.2. Model, simplifications and ansatz 
3.2.1. Neutral unpaired three-flavor quark matter 

The analysis of beta-equilibrated neutral three- flavor quark matter in an NJL model 
has been discussed in detail elsewhere (For example, Ref. EJ, and here we just sum- 
marize the main points. We assume that the up and down quarks are massless. The 
strange quark mass is a parameter M s . We couple a chemical potential (i to quark 
number, and always work at \i = 500 MeV. To ensure neutrality under color and 
electromagnctism we couple a chemical (electrostatic) potential \i e to negative elec- 
tric charge, and chemical (color-electrostatic) potentials fi^ and /ig to the diagonal 
generators of SU (3) co i OT . In full QCD these gauged chemical potentials would be the 
zcroth components of the corres ponding gauge fields, and would naturally be forced 
to their neutrality values ESlZZI i n the NJL context, we enforce neutrality explicitly 
by choosing them to satisfy the neutrality conditions J^- = = = 0. The 
assumption of beta-equilibrium is built into the calculation via the fact that the 
only flavor-dependent chemical potential we introduce is /Lt e , ensuring for example 
that the chemical potentials of d and s quarks with the same color must be equal. 

The presence of a nonzero strange quark mass, combined with the requirements 
of equilibrium under weak interactions and charge neutrality, drives down the num- 
ber of strange quarks, and introduces electrons and additional down quarks. We 
assume that this is a small effect, i.e. we expand in powers of jLt e //x and M^/fi . 
Working to linear order, we can simply treat the strange quark mass as a shift of 
—Mg/(2fj,) in the strange quark chemical potential. (The corrections from higher- 
order terms in M 2 / /i 2 in an NJL analysis of a two-flavor crystalline color supercon- 
ductor have been evaluated and found to be smallP 

and we expect the same to 

be true here.) In this approximation we find the Fermi momenta of the quarks and 
electrons are given by 




, Me u 
M + y Pf = V 



2/ie 

3 



Pf = Me 
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and their free energy is 



3(^) 4 + 3( P ^) + ( P |Q 4 3 r** 2 / Aa 

^unpaired = ^ J P W * ~ M ~~ T 



To this order, electric neutrality requires 

M 2 

i— it- (7) 

and because we are working to lowest order in M 2 / /i 2 we need no longer be careful 
about the distinction between p^'s and /x's, as we can simply think of the three 
flavors of quarks as if they have chemical potentials 

/id = M« + 2 <5/Z3 
= Pf 

Ms = /•*« ^ 25/i 2 (8) 
where, using Q, = /i — ^ and 

M 2 

= Sfj,2 = — - = Sfi , (9) 
8/i 

with the choice of subscripts indicating that 2Sfj,2 is the splitting between the Fermi 
surfaces for quarks 1 and 3 and 25^3 is that between the Fermi surfaces for quarks 
1 and 2, identifying u, d, s with 1, 2, 3. 

3.2.2. BCS pairing and neutrality 

BCS pairing introduces qualitative changes into the analysis of neutrality! ^ 1 | 33 | 3 4|4| 
For example, in the CFL phase /z e = and /is is nonzero and of order M 2 / 1 \x. This 
arises because the construction of a phase in which BCS pairing occurs between 
fermions whose Fermi surface would be split in the absence of pairing can be de- 
scribed as follows. First, adjust the Fermi surfaces of those fermions that pair to 
make them equal. This costs a free energy price of order (5/x 2 /i 2 . And, it changes the 
relation between the chemical potentials and the particle numbers, meaning that the 
/x's required for neutrality can change qualitatively as happens in the CFL example. 
Second, pair. This yields a free energy benefit of order A 2 ,/! 2 , where Ao is the gap 
parameter describing the BCS pairing. Hence, BCS pairing will only occur if the 
attraction between the fermions is large enough that Ao > <5/i. In the CFL context, 
in which (ud), (us) and (ds) pairing is fighting against the splitting between the d, 
u and s Fermi surfaces described above, it turns out that CFL pairing can occur if 
Ao > 4:5/1 = A/ s 2 /(2/i)Pa criterion that is reduced somewhat by kaon condensation 
which acts to stabilize CFL pairingEHEl 

We are now considering quark matter at densities that are low enough (p, < 
A/ 2 /(2Aq)) that CFL pairing is not possible. The gap parameter Aq that would 
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characterize the CFL phase if and 8^ were zero is nevertheless an important 
scale in our problem, as it quantifies the strength of the attraction between quarks. 
Estimates of the magnitude of Ao are typically in the tens of MeV, perhaps as large 
as 100 MeV0We shall treat Ao as a parameter, and quote results for Ao = 25 MeV, 
although our results can easily be scaled to any value of Ao as long as the weak- 
coupling approximation Aq <C [i is respected!^ 



3.2.3. Crystalline color superconductivity in two-flavor quark matter 

Crystalline color superconductivity can be thought of as the answer to the question: 
"Is there a way to pair quarks at differing Fermi surfaces without first equalizing 
their Fermi momenta, given that doing so exacts a cost?" The answer is "Yes, but it 
requires Cooper pairs with nonzero total momentum." Ordinary BCS pairing pairs 
quarks with momenta p and — p, meaning that if the Fermi surfaces are split at 
most one member of a pair can be at its Fermi surface. In the crystalline color 
superconducting phase, pairs with total momentum 2q condense, meaning that one 
member of the p air has momentum p + q and the other has momentum — p + q 
for some p! ' ' Suppose for a moment that only u and d quarks pair, making the 
analyses of a two-flavor model found in Refs. dl EU E21 ESI EH EHl EHI EZI EH1 (and 
really going back to Ref. 172(1 valid. We sketch the results of this analysis in this 
subsection. 

The simplest "crystalline" phase is one in which only pairs with a single q 
condense, yielding a condensate 

(ip u (x)Cj 5 ip d (x)) oc Aexp(2«q • r) (10) 

that is modulated in space like a plane wave. (Here and throughout, we shall denote 
by r the spatial three- vector corresponding to the Lorentz four-vector x.) Assuming 
that A <C Sfj, <C /J, the energetically favored value of |q| = q t urns o ut to be q = r]5fi, 
where the proportionality constant n is given by n = 1.1997.E^i^! If r\ were 1, then 
the only choice of p for which a Cooper pair with momentum (— p + q, p + q) would 
describe two quarks each on their respective Fermi surfaces would correspond to a 
quark on the north pole of one Fermi surface and a quark on the south pole of the 
other. Instead, with v > 1, the quarks on each Fermi surface that can pair lie on one 
ring on each Fermi surface, the rings having opening angle 2cos~ 1 (l/r7) = 67.1°. 
The energetic calculation that determines 77 can be thought of as balancing the gain 
in pairing energy as r\ is increased beyond 1, allowing quarks on larger rings to 
pair, against the kinetic energy cost of Cooper pairs with greater total momentum. 
If the A/8fi — » Ginzburg-Landau limit is not assumed, the pairing rings change 
from circular lines on the Fermi surfaces into ribbons of thickness ~ A and angular- 
extent ~ A/Sfi. 

After solving a gap equation for A and then evaluating the free energy of the 
phase with condensate 1(10(1 . one finds that this simplest "crystalline" phase is fa- 
vored over two-flavor quark matter with either no pairing or BCS pairing only within 
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a narrow window 



0.707 A 2SC < Sn < 0.754 A 2S c , 



(11) 



where A 2 sc is the gap parameter for the two- flavor phase with 2SC (2-flavor, 2-color) 
BCS pairing found at Sfi = 0. At the upper boundary of this window, A — * and one 
finds a second order phase transition between the crystalline and unpaired phases. 
At the lower boundary, there is a first order transition between the crystalline and 
BCS paired phases. The crystalline phase persists in the weak coupling limit only 
if <5/z/A 2 sc is held fixed, within the window while the standard weak-coupling 
limit A2SC/V — ^ is taken. Looking ahead to our context, and recalling that in 
three-flavor quark matter <5/i = Mf/(8/x), we see that at high densities one finds 
the CFL phase (which is the three-flavor quark matter BCS phase) and in some 
window of lower densities one finds a crystalline phase. In the vicinity of the second 
order transition, where A — > and in particular where A/5/j, — > and, consequently 
given (|1 1|) . A/A2SC ->0a Ginzburg-Landau expansion of the free energy order by 
order in powers of A is controlled. 

The Ginzburg-Landau analysis can then be applied to more complicated crys- 
tal structures in which Cooper pairs with sever al d ifferent q's, all with the same 
length but pointing in different directions, ariseP^ This analysis indicates that a 
face-centered cubic structure constructed as the sum of eight plane waves with q's 
pointing at the corners of a cube is favored, but it does not permit a quantitative 
evaluation of A(<5/z). The Ginzburg-Landau expansion of the free energy has terms 
that are quartic and sextic in A whose coefficients are both large in magnitude 
and negative. To this order, Q is not bounded from below. This means that the 
Ginzburg-Landau analysis predicts a strong first order phase transition between 
the crystalline and unpaired phase, at some Sfi significantly larger than 0.754 A 2 sc, 
meaning that the crystalline phase occurs over a range of dfi that is much wider than 
()11J) . but it precludes the quantitative evaluation of the Sfj, at which the transition 
occurs, of A, or of f2. 

In three-flavor quark matter, all the crystalline phases analyzed in Ref . 1521 have 
Ginzburg-Landau free energies with positive sextic coefficient, meaning that they 
can be used to evaluate A, fl and the location of the transition from unpaired quark 
matter to the crystalline phase with a postulated crystal structure. For the most 
favored crystal structures, we shall see that the window in parameter space in which 
they occur is given by ijljl. which is in no sense narrow. 

3.2.4. Crystalline color superconductivity in neutral three-flavor quark matter 

We now turn to three-flavor crystalline color superconductivity. We shall make 
weak coupling (namely Ao , <5/i <C /i) and Ginzburg-Landau (namely A -c Ao , Sfi) 
approximations throughout. The analysis of neutrality in three-flavor quark matter 
in a crystalline color superconducting phase is very simple in the Ginzburg-Landau 
limit in which A <C 5/j,: because the construction of this phase does not involve 
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rearranging any Fermi momenta prior to pairing, and because the assumption A -c 
5/j, implies that the pairing does not significantly change any number densities, 
neutrality is achieved with the same chemical potentials fx e = / (4/i) and ^3 = 
/is = as in unpaired quark matter, and with Fermi momenta given in Eqs. JHJ and 
(jHJ) as in unpaired quark matter. 

We consider a condensate of the form 

3 

(ip ia {x)C'Y 5 ip j p{x)) oc^2 H A / e 2iq ?- r e /Q)3 e / .y , (12) 
'=1 q?e{q 7 } 

where a, (3 are color indices, i, j are flavor indices, and where q", q% and q§ and A l5 
A2 and A3 are the wave vectors and gap parameters describing pairing between the 
(d, s), (u, s) and (it, d) quarks respectively, whose Fermi momenta are split by 2<5/ii, 
2<5/U2 and 2<5/t 3 respectively. From (jHJ, we see that S(j,% = Sfa = Sfj.i/2 = M^/{%fx). 
For each /, {q/} is a set of momentum vectors that define the periodic spatial 
modulation of the crystalline condensate describing pairing between the quarks 
whose flavor is not /, and whose color is not /. Our goal in this paper is to compare 
condensates with different choices of {qjj's, that is with different crystal structures. 
To shorten expressions, we will henceforth write XLa = yZ„a c! „ x . The condensate 
(|T2"|) has the color-flavor structure of the CFL condensate (obtained by setting all 
q's to zero) and is the natural generalization to nontrivial crystal structures of the 
condensate previously analyzed in Refs. l69M71l in which each {q 7 } contained only 
a single vector. 

In all our results, although not in the derivation of the Ginzburg-Landau ap- 
proximation itself in Ref . 1521 we shall make the further simplifying assumption that 
Ai = 0. Given that 5fi\ is twice <5/^2 or <5/i3, it seems reasonable that Ai <C A2, A3. 
We leave a quantitative investigation of condensates with Ai 7^ to future work. 



3.2.5. NJL Model and Mean-Field Approximation 

As discussed in Section 13.11 we shall work in a NJL model in which the quarks 
interact via a point-like four-quark interaction, with the quantum numbers of singlc- 
gluon exchange, analyzed in mean field theory. By this we mean that the interaction 
term is 

3 - 

^interaction = ~ gHlpT^ 1p) (ipT Au lp) , (13) 

The full expression for T Av is (T A ") ai j 3j = 7^ '(T A ) al3 5ij, where the T A are the color 
Gell-Mann matrices. The NJL coupling constant A has dimension -2, meaning that 
an ultraviolet cutoff A must be introduced as a second parameter in order to fully 
specify the interaction. Defining A as the restriction that momentum integrals be 
restricted to a shell around the Fermi surface, /i — A < |p| < /i + A, the CFL gap 
parameter can then be evaluated! ^®^ 

^2 



A = 2 3 A exp 



(14) 
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In the limit in which which A <C Ap, Su <^i n, all our results can be expressed 
in terms of Ao; neither A nor A appear!"*' This reflects the fact that in this limit 
the physics of interest is dominated by quarks near the Fermi surfaces, not near A, 
and so once Ao is used as the parameter describing the strength of the attraction 
between quarks. A is no longer visible; the cutoff A only appears in the relation 
between Ao and A, not in any comparison among different possible paired phases. 
In our numerical evaluations, we shall take fi = 500 MeV, A = 100 MeV, and adjust 
A to be such that A is 25 MeV. 

In the mean-field approximation, the interaction Lagrangian (|13|l takes the form 

^interaction = ^A(x)^ T + -1p T A(x)tp , (15) 

where A(x) is related to the diquark condensate by the relations 

A(x) = \\T A »{W T ){T Av ) T . (16) 
The ansatz <|12ll can now be made precise: we take 

A(a;) = A CF (x) ® (7 7 5 , (17) 

with 

3 

A CF (x) m , 03 =5^2 ^i^'^i^Hj • (18) 

1=1 qf 

Upon making the mean field approximation, the full Lagrangian is quadratic and can 
be written simply upon introducing two comp one nt Nambu-Gorkov spinors. Upon 
so doing, gap equations can easily be derived!^ Without further approximation, 
however, they are not tractable because without further approximation it is not 
consistent to choose finite sets {q/}. When several plane waves are present in the 
condensate, they induce an infinite tower of higher momentum condensates!^ The 
reason why the Ginzburg-Landau approximation, to which we now turn, is such a 
simplification is that it eliminates these higher harmonics. 



3.3. Ginzburg-Landau Approximation 

The form of the Ginzburg-Landau expansion of the free energy can be derived using 
only general arguments. 

We shall only consider crystal structures in which all the vectors qj in the 
crystal structure {q/} are "equivalent". By this we mean that a rigid rotation of 
the crystal structure can be found which maps any q^ to any other qj leaving the 
set {q/} invariant. (If we did not make this simplifying assumption, we would need 
to introduce different gap parameters A" for each vector in {q/}.) As explained 
in Section 13.2.41 the chemical potentials that maintain neutrality in three-flavor 
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crystalline color superconducting quark matter are the same as those in neutral 
unpaired three-flavor quark matter. Therefore, 

^crystalline = ^unpaired + f2(Ai, A2, A3) , (19) 

with ^unpaired given in © with Q, and with 0(0, 0, 0) = 0. Our task is to evaluate 
the condensation energy Q(Ai, A2, A3). Since our Lagrangian is baryon number 
conserving and contains no weak interactions, it is invariant under a global U(l) 
symmetry for each flavor. This means that f2 must be invariant under A/ — ► e^ 1 A/ 
for each I, meaning that each of the three A/'s can only appear in the combination 
AjAj. (Of course, the ground state can and does break these U(l) symmetries 
spontaneously; what we need in the argument we are making here is only that 
they are not explicitly broken in the Lagrangian.) We conclude that if we expand 
f2(Ai, A2, A3) in powers of the A/'s up to sextic order, it must take the form 

n({Aj}) = ^ E p ^ A *i A ' + \ ( E A(AJA/) 2 + E Pu a; a / a}a j ] 
ii \ 1 i>j ) 

+ 7/(A/A/) 3 + Y/TiJJ A *i A/A}A, 7 A}A, 7 + 7l23 AJ Ai A* A 2 A* A 3 ^J , 

(20) 

where we have made various notational choices for convenience. The overall prefac- 
tor of 2fi 2 /it 2 is the density of states at the Fermi surface of unpaired quark matter 
with M s = 0; it is convenient to define all the coefficients in the Ginzburg-Landau 
expansion of the free energy relative to this. We have defined Pj = dim{q/}, the 
number of plane waves in the crystal structure for the condensate describing pairing 
between quarks whose flavor and color are not /. Writing the prcfactor Pj multiply- 
ing the quadratic term and writing the factors of ^ and | multiplying the quartic 
and sextic terms ensures that the a/, /?/ and 77 coefficients are defined the same 
way as in Ref. EH The form of the Ginzburg-Landau expansion (|20fl is model inde- 
pendent, whereas the expressions for the coefficients a/, Pi, flu, 7/, 7/,/j and 7123 
for a given ansatz for the crystal structure are model-dependent. In Sections IV 
and V of Ref. [S3 these coefficients are evaluated by deriving the Ginzburg-Landau 
approximation to this model. 

We see in Eq. (|20fl that there are some coefficients — namely aj, j3i and 7/ 
— which multiply polynomials involving only a single Aj. Suppose that we keep a 
single A/ nonzero, setting the other two to zero. This reduces the problem to one 
with two-flavor pairing only, and the Ginzburg-Landau coefficients for this problem 
have been calculated for many different crystal structures in Ref. ESI We can then 
immediately use these coefficients, called a, (3 and 7 in Ref. 1651 to determine our 
a/, (3 1 and 7/. Using aj as an example, we conclude that 

ai = a( qiM = -1 + '£ bg ( P^-) - J log f -rr^H^) . (2D 



2qi \qi-8m) 2 \A{qf - 8$) 
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where (5/_ij is the splitting between the Fermi surfaces of the quarks with the two 
flavors other than / and qi = |qf| is the length of the q- vectors in the set {q/}. 
(We shall see momentarily why all have the same length.) In (|2I|) . A2SC is the gap 
parameter in the BCS state obtained with 8fij — and Aj nonzero with the other 
two gap parameters set to zero. Assuming that Ao <C /-t, this gap parameter for 2SC 
(2-fiavor, 2-color) BCS pairing is given b} ^M 



*2SC 



23A . (22) 



In the Ginzburg-Landau approximation, in which the Aj are assumed to be small, 
we must first minimize the quadratic contribution to the free energy, before proceed- 
ing to investigate the consequences of the quartic and sextic contributions. Notice 
that depends on the cutoff A and the NJL coupling constant A only through 
A2SC, and depends only on the ratios qi/Sfii and i5^//A2sc- <Xr is negative for 
Sfii /A 2 sc < 0.754, and for a give n value of this ratio for which ai < 0, ai is most 
negative and l|21|) is minimized for 



mmm 



qi = rytf/ij with 77 = 1.1997 , (23) 
where 77 is defined as the solution to 

1 log (2 + 1) = 1 . (24) 



2r? \t]-1 

We therefore set qi = rj S/ii henceforth and upon so doing find that 121(1 becomes 



Minimizing ai has fixed the length of all the q-vectors in the set {q/}, thus elimi- 
nating the possibility of higher harmonics. 

It is helpful to imagine the (three) sets {qj} as representing the vertices of 
(three) polyhedra in momentum space. By minimizing ai, we have learned that each 
polyhedron {q 7 } can be inscribed in a sphere of radius r/Sfij. From the quadratic 
contribution to the free energy, we do not learn anything about what shape poly- 
hedra are preferable. In fact, the quadratic analysis in isolation would indicate that 
if ai < (which happens for Sfii < 0.754 A2sc) then modes with arbitarily many 
different q/'s should condense. It is the quartic and sextic coefficients that describe 
the interaction among the modes, and hence control what shape polyhedra are in 
fact preferable. 

The quartic and sextic coefficients Pi and 7/ can also be taken directly from the 
two-flavor results of Ref. Note that these coefficients are independent of A2SC 
and the NJL coupling co nstant and cutoff, as long as the weak-coupling approxi- 
mation A 2 sc <C /i is valid Given this, and given the only dimensionful 
quantity on which they can depend is Sfii. They are therefore given by Pi = P/Sfij 
and 7/ = j/S/ij where Pi and 77 are dimensionless quantities depending only on 
the directions of the vectors in the set {q/}- P and 7 have been evaluated for many 
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crystal structures in Ref. ESI resulting in two qualitative conclusions. Recall that, 
as reviewed in Section ^.2. 31 the presence of a condensate with some qf corresponds 
to pairing on a ring on each Fermi surface with opening angle 67.1°. The first qual- 
itative conclusion is that any crystal structure in which there are two qj's whose 
pairing rings intersect has very large, positive, values of both and 7, meaning that 
it is strongly disfavored. The second conclusion is that regular structures, those in 
which there are many ways of adding four or six qj's to form closed figures in 
momentum space, are favored. Consequently, the favored crystal structure in the 
two- flavor case has 8 q" 's pointing towards the corners of a cubeP Choosing the 
polyhedron in momentum space to be a cube yields a face-centered cubic modulation 
of the condensate in position space. 

Because the 0i and 77 coefficients in our problem can be taken over directly 
from the two-flavor analysis, we can expect that it will be unfavorable for any of 
the three sets {q 7 } to have more than eight vectors, or to have any vectors closer 
together than 67.1°. It can be proved that 0jj is always positive, and in all cases 
investigated to date 777 j also turns out to be positive!^ This means that we know 
of no exceptions to the rule that if a particular {q 7 } is unfavorable as a two-flavor 
crystal structure, then any three-flavor condensate in which this set of q-vectors 
describes either the Ai, A 2 or A 3 crystal structure is also disfavored. 

The coefficients 0u and 777J cannot be read off from a two-flavor analysis 
because they multiply terms involving more than one A7 and hence describe the in- 
teraction between the three different A/'s. Before evaluating the expressions for the 
coefficients we make the further simpifying assumption that Ai = 0, because the 
separation S^i between the d and s Fermi surfaces is twice as large as that between 
either and the intervening u Fermi surface. This simplifies Ij20|l considerably, elimi- 
nating the 7123 term and all the 0u and 777 j terms except ^32, 7223 and 7332- We 
further simplify the problem by focussing on crystal structures for which {q2/ and 
{q?} are "exchange symmetric", meaning that there is a sequence of rigid rotations 
and reflections which when applied to all the vectors in {q 2 } and {q 3 } together 
has the effect of exchanging {^2} and {^3}- If we choose an exchange symmetric 
crystal structure, upon making the approximation that 8^2 = 8^3 and restricting 
our attention to solutions with A2 = A3 we have the further simplification that 

7322 = 7233- 

Upon making these simplifying assumptions, the only dimensionful quantity on 
which 032 and 7322 can depend is 6[i, meaning that ,$32 = ,$32/ '8[i 2 and 7322 = 
j322/5[i 4 y^ Here, 032 and 7322 are dimcnsionlcss numbers, depending only on the 
shape and relative orientation of the polyhedra {q2} and {qj}, which must be 
evaluated for each crystal structure. Doing so requires evaluating one loop Fcynman 
diagrams with 4 or 6 insertions of A/'s. Each insertion of A7 (AJ) adds (subtracts) 
momentum 2q" for some a, meaning that the calculation consists of a bookkeeping 
task (determining which combinations of 4 or 6 q°'s are allowed) that grows rapidly 
in complexity with the compexity of the crystal structure and a loop integration 
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that is nontrivial because the momentum in the propagator changes after each 
insertion. See Ref. 1521 where this calculation is carried out explicitly for 1 1 crystal 
structures in a mean-field NJL model upon making the weak coupling (Ao, Su <C (J-) 
approximation. Within this approximation, neither the NJL cutoff nor the NJL 
coupling constant appear in any quartic or higher Ginzburg-Landau coefficient, and 
they influence a only through A2SC ex Ao. Hence, the details of the model do not 
matter as long as one thinks of Aq as a parameter, kept <C u. 



3.4. General results 

Upon assuming that Ai = 0, assuming that the crystal structure is exchange sym- 
metric, and making the approximation that Su2 = Sfi3 = 5/.J, = M 2 /(8u), the free 
energy H20JI reduces to 



^(A 2 ,A 3 ) 



(/3(A 



1 1 , 



f 1 

f 7322 (A|Af 



P32MM 



(26) 



where (3, 7, P32 and 7322 are the dimensionless constants that can be calculated for 
each crystal structure,"**' and where the 6 //-dependence of a is given by Eq. I|25|l. 
It is easy to show that extrema of i7(A2,Aa) in (As, A.s)-space must either have 
A2 = A3 = A, or have one of A2 and A3 vanishing!"*' The latter class of extrema are 
two-flavor crystalline phases. We are interested in the solutions with A2 = A3 = A. 
(Furthermore, the three-flavor crystalline phases with A2 = A3 = A are electrically 
neutral whereas the two-flavor solutions in which only one of the A's is nonzero are 
not.ES) 

Setting A2 = A3 = A, the free energy becomes 



n(A) = 

where we have defined 



2Pa(5u)A 2 
Pes = W - 

Teff = 27 - 



2<V 

■ Pz2 
27322 



■P 



efi 



3(5/i' 



■7off 



(27) 



(28) 



We have arrived at a familiar-looking sextic order Ginzburg-Landau fre e en ergy 
function, whose coefficients can be evaluated for specific crystal structures!^ 

If Peg and 7eff are both positive, the free energy l|2Tjl describes a second order 
phase transition between the crystalline color superconducting phase and the normal 
phase at the Sfi at which a(5u) changes sign. From 1)25(1 . this critical point occurs 
where 5u = 0.754 A2S0 In plotting our results, we will take the CFL gap to be 
A = 25 McV, making A 2SC = 2 1 / 3 A = 31.5 MeV. Recalling that 5a = A/ s 2 /(8^), 
this puts the second order phase transition at 



Ml 
a 



6.03 A 2S c = 7.60 A = 190.0 MeV 



(29) 



a=0 
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(The authors of Refs. EH 1711 neglected to notice that it is A 2 sc: rather than the 
CFL gap Ao, that occurs in Eqs. 1(21(1 and 1(25(1 and therefore controls the Sfi at 
which a = 0. In analyzing the crystalline phase in isolation, this is immaterial 
since either Ao or A2SC could be taken as the parameter defining the strength of 
the interaction between quarks. However, we shall compare the free energies of the 
CFL, gCFL and crystalline phases, and in making this comparison it is important 
to take into account that A 2 sc = 2 1 / 3 A .) For values of M^/fi that are smaller 
than 1(29(1 (that is, lower densities), a < and the free energy is minimized by a 
nonzero A = A m j n and thus describes a crystalline color superconducting phase. 

If /3 c ff < and 7 c fi > 0, then the free energy 1(27(1 describes a first order phase 
transition between unpaired and crystalline quark matter occurring at a = a* = 
3 /3 2 ff /(32 P%g). At this positive value of a, the function f2(A) has a minimum at 
A = with O = 0, initially rises quadratically with increasing A, and is then turned 
back downward by the negative quartic term before being turned back upwards 
again by the positive sextic term, yielding a second minimum at 



also with f2 = 0, which describes a crystalline color superconducting phase. For 
a < a», the crystalline phase is favored over unpaired quark matter. Eq. ((25|> 
determines the value of Sfi, and hence M% / u, at which a — a* and the first order 
phase transition occurs. If a* <C 1, the transition occurs at a value of Mjy/x that 
is greater than 1)29(1 by a factor (1 + a*). Sec Fig. 0in Section l!T~71 for an explicit 
example of plots of f2 versus A for various values of a for one of the crystal structures 
that we describe. 

A necessary condition for the Ginzburg-Landau approximation to be quanti- 
tatively reliable is that the sextic term in the free energy is small in magnitude 
compared to the quartic, meaning that A 2 <C d/j, 2 \/3 e s/%g\. If the transition be- 
tween the unpaired and crystalline phases is second order, then this condition is 
satisfied close enough to the transition where A — > 0. However, if f3 e s < and 
%ff > 0, making the transition first order, we see from 1(30(1 that at the first or- 
der transition itself A is large enough to make the quantitative application of the 
Ginzburg-Landau approximation marginal. This is a familiar result, coming about 
whenever a Ginzburg-Landau approximation predicts a first order phase transition 
because at the first order phase transition the quartic and sextic terms are balanced 
against each other. We shall find that our Ginzburg-Landau analysis predicts a first 
order phase transition; knowing that it is therefore at the edge of its quantitative 
reliability, we shall focus on qualitative conclusions. 

3.5. Two plane wave structure 

We begin with the simplest three-flavor "crystal" structure in which {q 2 } and {q 3 } 
each contain only a single vector, yielding a condensate in which the (us) and 




(30) 
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(ud) condensates are each plane waves. This simple condensate yields a qualitative 
lesson which prov es e xtremely helpful in winnowing the space of multiple plane wave 
crystal structures!^ 

For this simple "crystal" structure, all t he coeffficients in the Ginzburg-Landau 
free energy can be evaluated analytically!^! The terms that occur in the three- flavor 
case but not in the two-flavor case, namely P32 and 7322, describe the interaction 
between the two condensates and hence depend on the angle <fr between q 2 and q 3 . 
For any angle </>, both ^32 and 7322 are positive, both increase monotonically with 
<j>, and both diverge as <fi — ► 7r!^ This tells us that within this two plane wave 
ansatz, the most favorable orientation is <j> = 0, namely q 2 || q 3 . Making this choice 
yields the smallest possible /3 e g and %s within this ansatz, and hence the largest 
possible A and condensation energy, again within this ansatz. The divergence at 
(ft — * 7r tells us that choosing q 2 and q 3 precisely antiparallel exacts an infinite free 
energy price in the combined Ginzburg-Landau and weak-coupling limit in which 
A <C Sfj,, Ao <C /it, meaning that in this limit if we chose <f> = tt we find A = 0. Away 
from the Ginzburg-Landau limit, when the pairing rings on the Fermi surfaces 
widen into bands, choosing <fi = it exacts a finite price meaning that A is nonzero 
but smaller than that for any other choice of </>. 

The high cost of choosing q 2 and q 3 precisely antiparallel can be understood 
qualitatively as arising from the fact that in this case the ring of states on the 
it-quark Fermi surface that "want to" pair with d-quarks coincides precisely with 
the ring that "wants to" pair with s-quarks!^ The simple two plane wave ansatz 
has been analyzed in the same NJL model that we employ upon making the weak- 
coupling approximation but without making a Ginzburg-Landau approximation!^ 
All the qualitative lessons that we have learned from the Ginzburg-Landau approx- 
imation remain valid, including the favorability of the choice cj> = 0, but we learn 
further that in th e tw o plane wave case the Ginzburg-Landau approximation always 
underestimates 

The analysis of the simple two plane wave "crystal" structure, together with the 
observation that in more complicated crystal structures with more than one vector 
in {q 2 } and {q 3 } the Ginzburg-Landau coefficient /3 32 (7322) is given in whole (in 
part) by a sum of many two plane wave contributions, yields an important lesson 
for how to constr uct favorable crystal structures for three-flavor crystalline color 
superconductivity!^ {q 2 } and {q 3 } should be rotated with respect to each other 
in a way that best keeps vectors in one set away from the antipodes of vectors in 
the other set. Crystal structures in which any of the vectors in {q 2 } are close to 
antiparallel to any of the vectors in {q 3 } are strongly disfavored because if a vector 
in {q 2 } is antiparallel (or close to antiparallel) to one in {q 3 }, this yields infinite 
(or merely large) positive contributions to P32 and to 7322 and hence to (3 e g and 
Teg , meaning that there is an infinite (or large) free energy penalty for A ^ 0. 

Summarizing, we have arrived at two rules for constructing favorable crystal 
structures for three-flavor crystalline color superconductivity. First, the sets {q 2 } 



February 2, 2008 6:19 WSPC/Trim Size: 9.75in x 6.5in for Review Volume beyond'BCS 




Fig. 5. The momenta that constitute the CubeX (left) and 2Cube45z (right) crystal structures. 
The set q3 of (ud) pairing momenta is given by the solid-tipped (red online) arrows. The set q2 
of (us) pairing momenta is given by the hollow- tipped (blue online) arrows. 

and {q 3 } should each be chosen to yield crystal structures which, seen as separate 
two- flavor crystalline phases, are as favorable as possible. In Section 13.2.31 we have 
reviewed how this should be done and the conclusion that the most favored {Qo| 
or {q 3 } in isolation consists of eight vectors pointing at the corners of a cube!^ 
Second, the new addition in the three-flavor case is the qualitative principle that 
{q 2 } and {q 3 } should be rotated with respect to each other in such a way as to 
best keep vectors in one set away from the antipodes of vectors in the other set. 



3.6. Multiple plane waves 

Rajagopal and Sharma have analyzed 11 different crystal structures in each case 
calculating f3 and 7 and ^32 and 7322, and hence /3 e ff and %s that specify the free 
energy as in (|28J) . The 11 structures allow one to make many pairwise comparisons 
that test the two qualitative principles described in Section 13.51 There are some 
instances of two structures which differ only in the relative orientation of {q 2 } and 
{q 3 } and in these cases the structure in which vectors in {q 2 } get closer to the 
antipodes of vectors in {q 3 } are always disfavored. And, there are some instances 
where the smallest angle between a vector in {q 2 } and the antipodes of a vector in 
{q 3 } are the same for two different crystal structures, and in these cases the one 
with the more favorable two-flavor structure is more favorable. These considerations, 
together with explicit calculations, indicate that two structures, which we denote 
"2Cube45z" and "CubeX" , are particularly favorable. 

The 2Cube45z crystal momenta are illustrated in Fig. [5] right panel; {q 2 } and 
{q 3 } each contain eight vectors pointing at the corners of a cube. If we orient 
{q 2 } so that its vectors point in the (±1,±1,±1) directions, then {q 3 } is rotated 
relative to {q 2 } by 45° about the z-axis. In this crystal structure, the (ud) and (us) 
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condensates are each face-centered cubi c, t he structure which in isolation is the 
most favored two-flavor crystal structure!"*' The relative rotation has been chosen 
to maximize the separation between any vector in {q 2 } and the nearest antipodes of 
a vector in {q 3 }. The color-flavor and position space dependence of the condensate, 
defined in l|17|) and (|18|1 . is given by 



AcF(x)ai,0j =£2aj3£2ij 2 A 



cos — [x + y + z) + cos — \—x + y + z) 
a a 



+ cos — (x — y + z ) + cos — ( — x — y + z) 
a a 



cos — {\Plx + z] + cos — (V%y + z 
a \ / a \ 



(31) 



+ cos — ^-V2y + z^j + cos — (—V2x + z^j 



where a and (3 (i and j) are color (flavor) indices and where 



V3tt _ 4.536 
q $[1 



1.764 M? 



(32) 



is the lattice spacing of the face-centered cubic crystal structure. For example, with 
A/ s 2 //i = 100,150,200 MeV the lattice spacing is a = 72,48,36 fin. Eq. (SJ can 
equivalently be written as Agf^cm,/^ = ^lafi^ ;A 2 (r) + e 3Q/ 3e 3 y A 3 (r), with (jSU 
providing the expressions for A 2 (r) and A 3 (r). A three-dimensional contour plot 
that can be seen as depicting either A 2 (r) or A 3 (r) separately can be found in 
Ref. ESI We have not found an informative way of depicting the entire condensate 
in a single contour plot. Note also that in (|31|l and below in our description of 
the CubeX phase, we make an arbitrary choice for the relative position of A 3 (r) 
and A 2 (r). We show in Ref. E21 that one can be translated relative to the other 
at no cost in free energy. Of course, rotating one relative to the other changes the 
Ginzburg-Landau coefficients /3 32 and 7 322 and hence the free energy. 

We now turn to a description of the CubeX crystal structure. We arrive at this 
structure by reducing the number of vectors in {q 2 } and {q 3 }. This worsens the 
two-flavor free energy of each condensate separately, but allows vectors in {q 2 } to 
be kept farther away from the antipodes of vectors in {q 3 }. We do not claim to 
have analyzed all structures obtainable in this way, but we have found one and only 
one which has a condensation energy comparable in magnitude to the 2Cube45z 
structure. The CubeX structure is illustrated in Fig. [SJ left panel; {q 2 } and {q 3 } 
each contain four vectors forming a rectangle. The eight vectors together point 
toward the corners of a cube. The two rectangles intersect to look like an "X" if 
viewed end-on. In the CubeX phase, the color-flavor and position space dependence 
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Fig. 6. The CubeX crystal structure of Eq. The fi gure extends from to a/2 in the 

x, y and z directions. Both A2(r) and A3(r) vanish at the horizontal plane. A2<r) vanishes 
on the darker vertical planes, and A3(r) vanishes on the lighter vertical planes. On the upper 
(lower) dark cylinders and the lower (upper) two small corners of dark cylinders, Aj(r) = +3.3A 
(Aa(r) = — 3.3A). On the upper (lower) lighter cylinders and the lower (upper) two small corners 
of lighter cylinders, Aa(r) = — 3.3A (As(r) = +3.3A). Note that the largest value of A/ (r) | is 
4A, occurring along lines at the centers of the cylinders. The lattice spacing is a when one takes 
into account the signs of the condensates; if one looks only at |A/(r)|, the lattice spacing is a/2, a 
is given in 1321 . In 1331 and hence in this figure, we have made a particular choice for the relative 
position of As(r) versus A2(r). In fact, one can be translated relative to the other with no cost in 
free energy. 



of the condensate is given by 



+ e3Q/3 £ 3y 2A 



2tt 2tt . . 

cos — (x + y + z) + cos — ( — x — y + z) 
a a 



cos — {—x + y + z) + cos — (x — y + z) 
a a 



(33) 



We provide a depiction of this condensate in Fig. 

We are confident that 2Cube45z is the most favorable structure obtained by 
rotating one cube relative to another. We are not as confident that CubeX is the best 
possible structure with fewer than 8+8 vectors. However, the two most favorable 
structures that we have found, 2Cube45z and CubeX, are impressively robust and 
make the case that three-flavor crystalline color superconducting phases are the 
ground state of cold quark matter over a wide range of densities. If even better 
crystal structures can be found, this will only further strengthen this case. 
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Fig. 7. Free energy Q vs. A for the CubeX crystal structure, described in Section 13.61 at four 
values of M^/fi. From top curve to bottom curve, as judged from the left half of the figure, the 
curves are M^/^i = 240, 218.61, 190, and 120 MeV, corresponding to a = 0.233, 0.140, 0, -0.460. 
The first order phase transition occurs at M% / fj, = 218.61 MeV. The values of A and C at the 
minima of curves like these arc what we plot in Figs. HI and |21 

3.7. Free energy comparisons and conclusions 

The gap parameter A and free energy 51(A) can be evaluated for all the crystal 
structures whose Ginzburg-Landau coefficients have been determined !^ For a given 
crystal structure, Q(A) is given by Eq. with (3 e g and 7 e ff taken from Table 

II of Ref. The quadratic coefficient a is related to 5 a by Eq. . Recall that 
we have made the approximation that <5/i2 = SfJ>3 = 5a = Mg/(8/x), valid up to 
corrections of order Mf//z 4 . At any value of M% / u, we can evaluate a(5a) and 
hence 51(A), determine A by minimizing Q, and finally evaluate the free energy £1 
at the minimum. In Fig. we give an example of fi(A) for various M^/a for one 
crystal structure with a first order phase transition (CubeX), illustrating how the 
first order phase transition is found, and how the A solving the gap equations — i.e. 
minimizing Q — is found. In Figs.nand[21 we plot A and f2 at the minimum versus 
Mg/u for the two most favorable crystal structures that we have found, namely the 
CubeX and 2Cube45z described in Section 

In Figs. U and O we have chosen the interaction strength between quarks such 
that the CFL gap parameter at M s = is Ao = 25 MeV. However, our results 
for both the gap parameters and the condensation energy for any of the crystalline 
phases can easily be scaled to any value of Ao. Recall that the quartic and sextic 
coefficients in the Ginzburg-Landau free energy do not depend on Ao. And, re- 
call from Eq. I|25|l that Ao enters a only through the combination A2sc/<^A t J where 
A 2S c = 23A and 5a = M s 2 /(8^). This means that if we pick a Ao ^ 25 MeV, 
the curves describing the gap parameters for the crystalline phases in Fig. ^ are 
precisely unchanged if we rescale both the vertical and horizontal axes propor- 
tional to Ao/25 MeV. In the case of Fig. |3 the vertical axis must be rescaled by 
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(Ao/25 McV) 2 . Of course, the weak-coupling approximation Ao <C /U, which we 
have used to simplify the calculation of the Ginzburg-Landau coefficients in Rcf . 
will break down if we scale Ao to be too large. We cannot evaluate up to what 
Ao we can scale our results reliably without doing a calculation that goes beyond 
the weak-coupling limit. However, such calculations have been done for the gCFL 
phase in Ref. 1441 where it turns out that the gaps and condensation energies plotted 
in Figs, n an d HI scale with Ao and A§ to good accuracy for Ao < 40 McV with 
/i = 500 MeV, but the scaling is significantly less accurate for Ao = 100 MeV. Of 
course, for Ao as large as 100 MeV, any quark matter in a compact star is likely 
to be in the CFL phase. Less symmetrically paired quark matter, which our results 
suggest is in a crystalline color superconducting phase, will occur in compact stars 
only if Ao is smaller, in the range where our results can be expected to scale well. 

Fig. H can be used to evaluate the validity of the Ginzburg-Landau approxima- 
tion. The simplest criterion is to compare the A's for the crystalline phases to the 
CFL gap parameter Ao- This is the correct criterion in the vicinity of the 2nd order 
phase transition point, where dfi = M s 2 /(8/j) r* Ao as in (|29|) . Well to the left, it is 
more appropriate to compare the A's for the crystalline phase to 5fi = /(8/x). By 
either criterion, the Ginzburg-Landau approximation is at the edge of its domain 
of validity, a result which we expected based on the general arguments of Section 
13.41 Therefore, although we expect that the qualitative lessons that we have learned 
about the favorability of crystalline phases in three-flavor quark matter are valid, 
and expect that the relative favorability of the 2Cubc45z and CubeX structures and 
the qualitative size of their A and condensation energy are trustworthy, we do not 
expect quantitative reliability of our results. There is therefore strong motivation 
to analyze crystalline color superconducting quark matter with these two crystal 
structures without making a Ginzburg-Landau approximation. It will be very inter- 
esting to see whether the Ginzburg-Landau approximation underestimates A and 
the condensation energy for the crystalline phases with CubeX and 2Cube45z crys- 
tal structures, as it does for the much simpler structure in which A2(r) and Aa(r) 
are each single plane waves !^ 

Fig. |3 makes it clear that the three-flavor crystalline color superconducting 
phases with the two most favored crystal structures that have been found are ro- 
bust by any measure. Their condensation energies reach about half that of the CFL 
phase at M s = 0. Correspondingly, these two crystal structures are favored over the 
wide range of j [i seen in Fig. [5] and given in Eq. . 

Taken literally, Fig. [2] indicates that within the regime (J3J of the phase diagram 
occupied by crystalline color superconducting quark matter, the 2Cube45z phase 
is favored at lower densities and the CubeX phase is favored at higher densities. 
Although we do have qualitative arguments why 2Cube45z and CubeX are favored 
over other phases, we have no qualitative argument why one should be favored over 
the other. Also, in this context the Ginzburg-Landau approximation is not reliable 
enough for us to conclude that one phase is favored at higher densities while the 
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other is favored at lower ones. The proper conclusion from these results is that 
2Cubc45z and CubeX are the two most favorable phases we have found, that both 
are robust, that the crystalline color superconducting phase of three- flavor quark 
matter with one crystal structure or the other occupies a wide swath of the QCD 
phase diagram, and that their free energies are similar enough to each other that 
a definitive comparison of their free energies will require a calculation that goes 
beyond the Ginzburg-Landau approximation that we have used here. 

3.8. Implications and Future Work 

The purpose of this calculation has been to find candidates for the true ground 
state of quark matter in the "gCFL window" , the range of values of / [i where 
calculations that are limited to isotropic phases predict that the gCFL phase is 
favored. Fig. [5] shows that, according to a Ginzburg-Landau theory keeping terms 
up to scxtic order, in most of this range the gCFL phase can be replaced by a 
much more favorable three-flavor crystalline color superconducting phase and that 
these crystalline phases may be favored over a wide range of M% / fj, above the gCFL 
window also. Fig. [21 also indicates that at the lower end of the range of M^/u, 
closest to the CFL— >gCFL transition, there is a narrow region where it is hard 
to find a crystalline color superconducting phase with lower free energy than the 
gCFL phase. This narrow region is thus the most likely place to find a ground state 
consisting of the gCFL phase augmented by current-carrying meson condensates 
described in Refs. 1141 1181 Except within this window, the crystalline color super- 
conducting phases with cither the CubeX or the 2Cube45z crystal structure provide 
an attractive resolution to the instability of the gCFL phase. 

Our Ginzburg-Landau calculation provides reason to speculate that most or 
possibly all of the quark matter core of a compact star might be in a crystalline 
color superconducting phase. For example, substituting the quite reasonable values 
A = 25 McV and M s = 250 MeV into Eq. (gj we find that the crystalline phases are 
favored over the range 240McV < u < 847MeV. Obviously at the lower density end 
of this range all quark matter phases are replaced by nuclear matter, and the high 
end of this range extends far beyond the values /j, ~ 500 MeV that are expected 
in the center of compact stars. We therefore find that with these choices for the 
parameters the crystalline phase covers the entire range of expected densities of 
stellar quark matter. Of course, if Ao is larger, say ~ 100 MeV, the entire quark 
matter core could be in the CFL phase. And, there are reasonable values of Ao 
and M s for which the outer layer of a possible quark matter core would be in a 
crystalline phase while the inner core would not. We do not know Ao and M s well 
enough to predict what phases of quark matter occur in compact stars. However, 
from our results it is clear that one of the possible configurations of a compact star, 
ultimately to be confirmed or refuted by astrophysical observations, would involve 
a quark matter core that is mostly or entirely in a crystalline color superconducting 
phase. 
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Now that we have two well-motivated ansatze for the crystal structure of three- 
flavor crystalline color superconducting quark matter and a good qualitative guide 
to the scale of A and fi(A), we should be able to make progress toward the cal- 
culation of astrophysically observable preperties of compact stars. For example, we 
would like to predict the effects of a crystalline color superconducting quark core 
on the rate at which neutron stars cool by neutrino emission and on the occurrence 
and phenomenology of pulsar glitches. 

For cooling, the relevant microphysical quantities are the specific heat and the 
neutrino emissivity. The specific heat of the crystalline phases will not be dramati- 
cally different from that of unpaired quark matter: it rises linearly with T because 
of the presence of gapless quark excitations at th e bo undaries of the regions in 
momentum space where there are unpaired quarks (In contrast, in the gCFL 
case the heat capacity is parametrically enhanced.^) However, the neutrino emis- 
sivity should turn out to be significantly suppressed relative to that in unpaired 
quark matter. The evaluation of the phase space for direct Urea neutrino emission 
from the CubeX and 2Cube45z phases will be a nontrivial calculation, given that 
thermally excited gapless quarks occur only on patches of the Fermi surfaces, sep- 
arated by the (many) pairing rings. (The direct Urea processes u + e — * s + v and 
s — > u + e + v require s, u and e to all be within T of a place in momentum space 
where they are gapless and at the same time to have p u + pe = p s to within T. 
Here, T ~ keV is very small compared to all the scales relevant to the description 
of the crystalline phase itself.) 

The idea that a region of a star made of quark matter in a crystalline phase 
could be the place where (some) glitches originate goes back to Ref. 1191 where it was 
suggested that a crystalline color superconductor condensate could pin rotational 
vortices. This leads to the possibility that the presence of crystalline quark matter 
in compact stars could be strongly constrained by comparing their predicted glitch 
phenomenology with what is actually observed. 

There are two key microphysical properties of crystalline quark matter that must 
be estimated before glitch phenomenology can be addressed: the pinning force on 
rotational vortices and the shear modulus of the crystal. In order for glitches to 
occur it is necessary for rotational vortices in some region of the star to be pinned 
and immobile while the spinning pulsar's angular velocity slows over years. But 
this can only occur if the pinning force is strong enough to lock the vortices to the 
crystal structure and if the structure resists deformation, i.e. has a large enough 
shear modulus. (The glitches themselves are then triggered by the catastrophic 
unpinning and motion of long-immobile vortices.) 

To estimate the pinning force we must analyze how the CubeX and 2Cubc45z 
phases respond when rotated. We expect vortices to form and be pinned at the 
intersections of the nodal planes at which condensates vanish. This analysis faces 
two complications. Firstly, baryon number current can be carried by gradients in 
the phase of either the (us) crystalline condensate or the (ud) condensate or both, 
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so we must determine which types of vortices are formed. Secondly, the vortex core 
size, 1/A, is only a factor of three to four smaller than the lattice spacing a. This 
means that the vortices cannot be thought of as pinned by an unchanged crystal; 
the vortices themselves will qualitatively deform the crystalline condensate in their 
vicinity. 

The second microphysical quantity that is required is the shear modulus of the 
crystal. This can be related to the coefficients in the lo w energy effective theory 
that describes the phonon modes of the crystal This effective theory has 
been analyzed, and its coefficients calculated, for the two-flavor crystalline color 
superconductor with face-centered cubic symmetry!^! Analyzing the phonons in 
the three-flavor crystalline color superconducting phases with the 2Cube45z and 
CubeX crystal structures is thus a priority for future work. 

The instability of the gCFL phase was a significant bend in the road to un- 
derstanding the properties of dense matter, but one that is now receding behind 
us. With the emergence of two robust crystalline phases, possibly accounting for 
as much as all of the quark matter that may occur, we can glimpse on the road 
ahead the microphysical calculations prerequisite to identifying observations that 
can be used to rule out (or in) the presence of quark matter in a crystalline color 
superconducting phase within neutron stars. 

4. Coda 

The project of delineating a plausible phase diagram for high-density quark mat- 
ter is still not complete. We have discussed some ideas for the "non-CFL" region, 
but th ere are others such as a suggested gluon condensation in two-flavor quark 
matterpHand deformation of t he Fermi surfaces (discussed so far only in non-beta- 
cquilibrated nuclear matteJ22l). It is very interesting to note that the problem of 
how a system with pairing responds to a stress that separates the chemical po- 
tentials of the pairing species is a very generic one, arising in condensed matter 
systems a nd c old atom systems as well as in quark mat ter. Recent work by Son and 
StephanoJ^l and Mannarelli, Nardulli and RuggierP^l on a two-species model char- 
acterized by a diluteness parameter and a splitting potential shows that between 
the BCS-paired region and the unpaired region in the phase diagram one should 
expect a translationally-broken region. In QCD this could correspond to a p-wave 
meson condensate or a crystalline color superconducting state. What is particularly 
exciting is that the technology of cold atom traps has advanced to the point where 
fermion superfluidity can now be seen in conditions where many of the important 
parameters can be manipulated, and it may soon be possible to investigate the re- 
sponse of the pairing to external stress under controlled experimental conditions. 
Pairing in a gas of ultracold fcrmionic atoms can never yield a complete analogue of 
pairing in dense quark matter, because the atoms are neutral and thus the require- 
ment of neutrality imposes no constraints, whereas we have seen that in the quark 
matter context imposing both color and electric neutrality is a crucial qualitative 
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driver for the physics, for example precluding phase separation which can easily 
happen in the cold atom system. However, even without yielding a complete ana- 
logue these experiments could answer qualitative questions of interest. For example, 
should they discover a translationally broken regime in the absence of rotation, it 
will be very interesting first to see what crystal structure emerges and second to 
study the response of such a phase to rotation. 
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